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Abstract 

We investigated the possibility of testing factorization hypothesis in non- 
leptonic exclusive decays of 5-meson. In particular, we considered the non- 
factorizable — > modes and —>■ (vr^, p^) known as 

well-factorizable modes. By taking the ratios D^*^^Di*^^)/B{B^ 

(| 7]-~, we found that under the present theoretical and experimental 
uncertainties there's no evidence for the breakdown of factorization description 
to heavy-heavy decays of the B meson. 
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I. INTRODUCTION 



Non-leptonic decays of heavy mesons are very important weak processes for the determination 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements and the understanding of the 
CP violation mechanism. The analysis of decay processes of S-mesons, now produced in large 
numbers in the 5-factories worldwide, will contribute to the stringent test of the Standard 
Model (SM) and search for new physics. Non-leptonic decays of S-mesons, however, are 
complicated processes due to inherent hadronic nature and final state interactions. 

A simple formulation of the decay amplitude, so called the naive factorization scheme 
[0,^, has been widely used without full theoretical justification. And its phenomenologi- 
cal extension, the generalized factorization scheme, with process-dependent quantities from 
penguin effects and non-factorizable contributions has been also widely used in the literature 
In this latter scheme, the non-factorizable effects are contained in the effective color 
number, A*"^^, which is a free parameter of the scheme; the value of A^^^ was adjusted 
to oo for D decays, and to 2 or 5 depending on the chiral structure of B decays. However, 
the A'^^'^-dependence of the Wilson coefficient in the effective Hamiltonian is different for 
each individual coefficient. For example, the coefficients related to color-favored processes, 
ai, 04, ag, ag are stable against the variation of A'^''', while those related to color-suppressed 
processes strongly depend on A"^''^ 0^]- This observation indicates that the non-factorizable 
contributions in the color-favored decays are negligible compared with the factorizable con- 
tributions. 

Recently much progress has been achieved towards understanding non-leptonic decay 
processes by separating out short distance physics from long distance effects in the well- 
defined manner; Beneke et al. proved the validity of factorization for the i?-meson decay 
amplitude in the context of perturbative QCD formalism. They showed that when a i?-meson 
decays weakly to a heavy meson and emits a light meson, the decay amplitude factorizes 
in the same form as the naive factorization formula, but with calculable coefficients in the 
heavy quark limit. In the case of 5-meson decaying to light meson rather than heavy 
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one, according to their formalizm, a contribution by hard spectator quark is added to the 
amphtude, therefore the total amphtude is still factorizable. However, for 5-decays to a 
heavy or light meson with emitting a heavy meson, their amplitudes are not written as 
factorized forms, since the color transparency arguments cannot be applied for such decays. 

Though the factorization of the decay amplitude for a 5-meson decaying to heavy-heavy 
mesons has not been justified, there have been many calculations using the factorized for- 
mula in the literature P,P,pil|. Within the naive factorization scheme, Luo and Rosner [|I^] 
calculated the branching ratios of the 5-meson decays, Z}(*)+D^*)~, after extracting 

the values of \Vch\ and the slope of the universal Isgur-Wise form factor p^, by comparing 
the decay rates of D^*^^{n^ , p^) with a differential distribution of 5^ D'^*^^l^ui 

measured by the CLEO Collaboration |Tl[|. They found theoretical predictions of the naive 
factorization approach are acceptable within present experimental errors. 

Here we test the generalized factorization scheme for the color-favored i?-meson decay 
to heavy-heavy mesons by comparing with the 5-decay to heavy-light mesons. The selected 
decay processes are B^ £)(*)+£){*)- for heavy-heavy and B^ D^*^^{n^ , p^) for heavy- 
light, whose experimental branching ratios are well known. Compared to the work of Luo 
and Rosner, in which the authors used the naive factorization scheme neglecting penguin 
effects, we here include penguin effects and take ratios of the decay rates to reduce the form 
factor dependence and cancel the CKM matrix elements. In the next section, we will give 
theoretical descriptions of B^ D^*^^{ti^ , p~) and B^ /}(*)+/}(*)- within generalized 
factorization scheme, and present numerical analyses. Section IH discusses experimental 
feasibility of our analyses and closes with a brief summary. 

II. THEORETICAL DESCRIPTION OF D(*)+(7r ,p ) AND B° Z)(*)+i:)^*^ 

WITHIN GENERALIZED FACTORIZATION APPROACH 

As previously mentioned, there has been general consensus on the applicability of factor- 
ization approach to the color-favored heavy-light decays B D^*\7i,p), and this has 
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been recently justified within perturbative QCD formalism 0. However, questions still 
remain about the applicability of factorization to heavy-heavy decay of 5-meson, such as 
B — i> D'^*^D^*\ We here investigate the validity of factorization hypothesis by taking ratios of 
branching fractions of presumably non-factorizable 5^ —>■ Z}(*)+/}(*)^ modes to factorizable 
# ^ L)W+(7r-,p-) modes. 

Based on the generalized factorization formalism, the decay amplitudes of our interest 
are expressed as 

A(b' ^ dW+M-) = ^\4,\/;,S(DWM)(M-|gV(l - l,)q\0){D^*^^\c^,{l - ^,)b\B'), (1) 

where q{q') = u{d) for M = tt, p and q{q') = c(s) for D^*\ The coefficient a includes penguin 
effects and possible non-factorizable contributions in the generalized factorization scheme. 
They are given, neglecting W-exchange diagram and using VtbV^* = —VcbV*^, as 



a{DDs) = ai 1 + " + 2 



04 + aio , ^06 + as 



ai ai {rrib - rric) {rric + m J 



KZ^WD:) = aifl + ^i±^V (2) 



ai 

where a^-'s represent conventional effective parameters defined as a2i = cf + cl^_i/Nf and 
a2i_i = c'^-i+c'^ / . Using the numerical values of a/s of Ref. 0, the effective parameters 
d defined above are related to ai by 

\d{B DDs)\ = 0.847ai, 
\d{B D*Ds)\ = 1.037ai, 
\a(B D^*^D*)\ = 0.962ai, (3) 

where the values are obtained by choosing N^^ = 2 for {V ~ A){V — A) interactions {i.e. for 
operators Oi^2,3,4,9,io) and A''^^ = 5 for (V — A)(y + A) interactions {i.e. for operators Os^gj.s) 
0. We note that the ratios, |a/ai|, are numerically very stable over different N^^ values; 
for example, the numerical deviations are less than a few % for N^^ = 2, 3, 5 and oo. From 



the relations in Eq. one can see that, at the amphtude level, the penguin contributions 
to fi^ D*^Dj decay (~ 3.7%) are much smaller than those for D^Dj mode 

(~ 15.3%). In fact the penguin effects on b'^ — D^D^ decay are not small enough to 
be simply neglected. As previously mentioned, the penguin effects are neglected in the 
analyses of Ref . |TD[ . We will show that the inclusion of the penguin effect in B D+D' 
mode improves substantially the theoretical prediction to the experimental value. For fi^ —* 
D*^Dj decay mode, the penguin contribution can be neglected. This difference of penguin 
contributions between the similar decay modes 5^ D^Dj and B^ D*^Dj is due to 
the different chiral structure of the final states: B —>■ D* transitions occur through axial 
vector currents, while B D through vector currents. 
The ratios 
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are given as 
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.fp J \Pc' J \ mp J \Ai 
where is the cm. momentum of the decay particles and we used |V"^*/Kjrf| = 1. Here 
the form factors have the following parameterization : 



{P'{P)K\P{P)) = ipp+p', - "^^^q^ Fi(g^) + "^^^-^qMq^ 



{V{p\e)\V,\P{p)) = ^ 6^.„;3e>W(g2) 

mp + my 



{V{p\e)\A,\P{p)) = i 



[mp + my)e^Ai(g2) "-^—{p + p')/.^2(g') 

mp-\-mv 



(10) 



where q = p - p' , Fi(0) = Fo(0), ^3(0) = Ao(0), 



and P, V denote the pseudoscalar and vector mesons, respectively. For B ^ V1V2 decay 
(see Eq. (|)), three form factors Ai{q^), ^2(9^), and V{q^) contribute. Here we factored out 
the dominant one Ai{q^) and the other two are put in the function H{q^) defined as 

H{q^) = {a-hxf + 2{l + cY), (11) 

with 



— mi — m2 2m^p^ 2mBPc 



2mim2 ' mim2{mB + TTLiY'' (m^ + mi''^' 



(12) 



A^^q^y ^ 

where mi (m2) is the mass of the vector meson Vi (V2). Using the above ratios Eqs.(|^)- 
(^), one can, in principle, test the validity of factorization without having dependence on 
CKM matrix elements. However, the analysis depends strongly on nonperturbative hadronic 
factors such as decay constants and form factors. B D^*^ transition form factors are 
rather well-constrained and the uncertainty in their ratios would be rather moderate. In 
the following numerical analysis, we consider three models for the form factors of i? ^ D'^*^ 
transitions: the Bauer-Stech-Wirbel (BSW) model 0, Melikhov/Stech ||12|, and relativistic 
light-front (LF) quark model ||13|| . 

Another uncertainty comes from decay constants, especially /^(*), which has presently 



large uncertainty. Particle Data Group report gives two distinct values depending on its 
decay modes: 



5 



= 194 ± 35 ± 20 ± 14 MeV from Ds fii^f,, 
/^+ = 309 ± 58 ± 33 ± 38 MeV from tu^. 



(14) 
(15) 



Recently, a rather interesting value appeared in Ref . [|r^ : 



= 323 ± 44 ± 12 ± 34 MeV from ^z/^, 



(16) 



which is obtained by measuring the branching fraction of Dg fiu^ relative to the branching 
fraction Dg ipn K^K^tt. We use the statistical average of the above three, 



/d, = 252 ± 31 MeV. 



(17) 



Then we get the theoretical predictions 



T^Ds/n = [3.38 ± 0.841] 
Ud'J^ = [0.92 ± 0.236] 
Uojn = [2.17 ±0.654] 
TZdi/p = [2.15 ±0.545] 



a{D7T) 
~a{DDl 



f. 



0.686 



0.252/ V 0.74 ; \F^^{ml) 
\ i'l ("^D^ 




:i^ 



where the quoted errors are based on our estimates of uncertainties in the form-factor model- 
dependence and in the decay constants /^{.). Here we assumed fa* = f^s for simplicity and 
used = 131 MeV and fp = 209 MeV [0]. As the ratios of a's are factored out, the numerical 
predictions of Eqs. ([T8|) correspond to those in the naive factorization approximation. As is 
shown, the main uncertainty comes from our ignorance on the decay constant • Within 
the generalized factorization scheme and by including penguin effects, the central values of 
the ratios are shifted to 



7^gf/^ = 2.43 ±0.61, 
n%^jp = 0.85 ± 0.22, 
^gf/^ = 2.33 ± 0.70, 
:^GF^^ = 2.00 ± 0.50, 



(19) 



where we used the exphcit numerical values for a of Eq. (|^). Considering the current exper- 
imental branching ratios for each decay mode [jlO|,|l^, one gets the following ratios: 



7^^P^ = 2.67 ± 1.061, 
7^^f/^ = 1.27 ±0.671, 
:^^P^ = 3.58 ± 1.138, 

= 2.16 ±0.817. (20) 

Comparing the ratios ([T8|) , ([ToD and (^), all the theoretical predictions are well within 
the present experimental constraints. We note that the inclusion of the penquin effects for 
5^ — > DjD^, which add a sizable contribution, improves the central value so that it is 
much closer to the experimental value. Although presently the experimental errors are too 
large to say anything definite, our analysis indicates that factorization hypothesis is still a 
good method for describing the 5-meson decaying to heavy-heavy mesons. Furthermore, 
one could even consider a possibility that the factorization may not be a consequence of 
only perturbative QCD, in contrast to the arguments of Ref. [0. Similar arguments are 
given in Ref. [|1^], in which the authors considered B D^*^X decays and expected non- 
factorization effects would grow with the invariant mass of the multi hadronic state X if 
the factorization is a consequence of perturbative QCD, but they found no such dependence 
on nix- 



III. DISCUSSIONS ON EXPERIMENTAL FEASIBILITY AND SUMMARY 

The comparison of (|19D and (^) will give a test of generalized factorization model that we 
considered in this paper. As it stands now, the two sets of values are consistent well within 
uncertainties. On the theoretical side, the biggest uncertainty is in the determination of 
meson decay constants /^c*), while on the experimental side the statistical errors of B{B — > 
give the largest uncertainty. Therefore, we need to improve the precision of such 
experimental measurements, for the method described in this paper to have any significance. 
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As of this writing, the combined data sample of BABAR and Belle experiments is 
more than 30 fb^^ [|l^. The peak instantaneous luminosity of each experiment is over 
3 X lO'^'^cm^^s^^, which corresponds to more than 1 fb^^/week for each. At this rate, we 
expect to have more than 100 fb~^ of data accumulated within a year. We will take this as 
our basis for considering experimental feasibility. 

Currently, the most precise measurement of /^j^ is obtained by the CLEO collabora- 
tion |T^ in Ds jJiu decays. Adding the errors in quadrature, they obtained /d^ = 280 ±48; 
17 % total uncertainty (7 % statistical) in /d^, corresponding to 34 % error in our calcula- 
tion of the ratios. With 100 fb~^ data sample from BABAR and Belle, which is more than 
20 times that of the existing result , the statistical error will be reduced to 1/ v^20 of 



Ref. |]T8|. The systematic errors may not go down as fast, but better understanding of every 
other aspect of the analysis will help reduce the systematic uncertainties. Assuming that 



the systematic error can be reduced to 1/3 of Ref. |T8|, the /^^ value will be determined 
to 5 % accuracy, hence resulting in 10 % error in our ratio. As for the form factor errors 
in the theoretical calculations, we hope that in the near future precision measurements in 
heavy-flavor physics processes from B or charm factories should help test and confirm the 
reliability of lattice QCD technique greatly. Then we may have much improved form factor 
errors. 

In the experimental measurements of branching ratios, B D^*^7i modes are measured 
with much better precision than B D^*^p modes. Similarly, D^*^Ds modes are determined 
with significantly higher precision than D^*^D* modes. Therefore, we expect that Rd^/tt 
and Rds/tt will be determined with higher precision than other ratios. The most recent and 
precise measurement of B{B D*+7r") = (2.81 ± 0.24 ± 0.05) x 10~^ is obtained by the 
CLEO collaboration with 3.1 fb~^ of data sample. With a data sample of 100 fb~^, we can 
measure this branching ratio with a precision at a few % level. On the other hand, the 
most precise measurement of B{B D*^Dj) also made by CLEO with 2.0 fb~^ of data, is 
(0.90 ±0.27 ±0.22) x 10"^. Again, with 100 fb"^ of data, the statistical error will be reduced 
to ~ 1/7 of its current value. If experimental systematic error can be reduced to ~ 40 % 
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of its current value, B{B D*^D^) will be determined to 10 % level. Therefore, it is 
likely that the experimental value of Rds/tt can be determined at the level of 10 % precision. 
Similar case can be made with Rd^/tt- 

Comparing (|T9|) and (pOf), we note that the experimental and theoretical values of Rds/tt 
show the biggest difference if we accept their central values. We also note that if both values 
can be determined within 10 % accuracy and if we assume that their central values stand 
as they are, then we will be able to see 3 a difference in Rd^/tt- In conclusion, we will have 
a good opportunity to test the generalized factorization scheme as discussed in this paper 
once we have ~ 100 fb~^ of data from the 5-factory experiments. 

To summarize, we have investigated the possibility of testing factorization hypoth- 
esis from non-leptonic exclusive decays of B meson into two meson final states. In 
particular, we considered the presumably non-factorizable 5^ /)(*)+/)(*)- modes and 
— s> D^*''^ (n^ , p^) known as well-factorizable modes. By taking the ratios B(B^ — > 
/)(*)+£)(*)- )/^(^^*^ ^ D*^*''+(7r^, p^)), the dependence on CKM matrix elements vanishes 
and some model-dependence on hadronic form-factors is reduced. We found that under 
the present theoretical and experimental uncertainties there's no evidence for breakdown of 
factorization description to heavy-heavy decays of the 5-meson. 
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